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Abstract 

Persistent current generation by small polarons embedded in a spin-density wave background 
in a two dimensional lattice is theoretically studied. When the embedded small polarons become 
cores of merons and antimerons (vortices in spin configuration with winding numbers +1 and — 1, 
respectively), a spin Berry phase arises from the spiral spin texture of the merons and antimerons, 
and generates a spin current around the polarons. Electric current is, however, absent due to 
the exact cancellation of the up- and down-spin electron currents. By addition of an imaginary 
magnetic flux that negates the fictitious magnetic flux from the spin Berry phase, the state with 
merons and antimerons is stabilized and non-zero electric current appears; this imaginary flux 
becomes a real one when side-by-side bipolaronic cores of merons and antimerons are formed, 
where the real flux arises spontaneously from the wave function for the motion of the system as 
a whole; thus, the bipolarons serve as pumps for the persistent electric current. It is also shown 
that a Nernst signal stems from thermal diffusion of the bipolarons. Implications of the present 
results to anomalous properties of cuprate superconductors, including the enhanced Nernst signal, 
Drude peak, and Fermi arc are also discussed. The present result suggests the appearance of the 
phase variable or "Nambu-Goldstone" mode for the superconductivity without the breakdown of 
the global U(l) gauge invariance in the cuprates. 

PACS numbers: 72.15.Gd,03.65.Ca,74.72.-h 
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I. INTRODUCTION 



F. Bloch once argued that a persistent current might arise from the absence of the limit, 

lim ^, (1) 

where AE is an increase of the total energy when a momentum in one direction P is in- 
creased by AP.— The momentum change AP is brought about by multiplying a function 
cxp (iAPXh' 1 ) on the wave function \1/ as 

■q,' = q e 1 ^, (2) 

where X denotes the coordinate of the center of gravity in the direction of P. Actually, the 
multiplied function is a wave function for the motion of the system as a whole. 
In classical mechanics, the limit corresponding to Eq. ([TJ is given by 

which is zero at the minimum of the Hamiltonian H; thus, a persistent current, X ^ 0, at 
the minimum is impossible. 

According to Bloch, however, in quantum mechanics, the limiting process in Eq. may 
not be permissible because of the discreteness of the momentum. Let us consider a system 
of N electrons confined in a narrow loop of circumference L, and take the x coordinate along 
the loop. The single-valuedness requirement on exp (iAPXh' 1 ) after the increase of x by L 
demands that AP must be 

AP = ^n, (4) 

where n is an integer. If L is macroscopically large, AP is practically continuous, thus, a 
result similar to the classical one will be obtained. However, if L is a microscopic length 
(for example, several atomic distance in a crystal), the discreteness is not negligible; then a 
result different form the classical one may arise. 

The above argument suggests that a persistent current may arise in a system with cavity 
regions of microscopic sizes. In the present work, we deal with such a system, namely, a 
system composed of conduction electrons and small polarons with the condition that the 
hopping motion of the latter is much slower than that of the former, thus, the latter provide 



lattice spacing size cavities for the motion of the former. Then, a macroscopic persistent 
current may be generated as a collection of circular currents around the small polarons. 

Actually, for the persistent current generation, the existence of microscopic cavities is 
not enough. A condition that stabilizes a current carrying state over a currentless state is 
necessary. We will show that such a condition is provided if the conduction electrons form 
a background spin density wave (SDW). Thus, the system considered in this work is the 
inhomogeneous system composed of the conduction electrons that form an SDW background 
and the small polarons that provide microscopic cavities. The inhomogeneity considered in 
this work may be relevant to cuprate superconductors, where the parent compounds of 
the cuprate superconductors are half-filled insulators with an antiferromagnetic SDW-like 
correlation; and doped holes are expected to form small polarons due to strong Coulomb 
interaction and strong vibronic interaction.— i 3 i 4 i 5 i 6 

When holes are doped in the cuprates, they are argued to become cores of merons (vortices 
in the spin configuration with winding number +1) or antimerons (vortices in the spin 
configuration with winding number — 1)^ and the existence of such objects explains the 
peak splitting observed in the neutron scattering experiment very naturally. The formation 
of merons and antimerons around the small polarons is the key phenomenon that leads 
to the appearance of a persistent current, which will be shown by employing a mean field 
approximation that leads to an SDW state in a Hubbard model. We will show that a spin 
current appears when the polarons become cores of merons and antimerons due to a spin 
Berry phase from the spiral spin texture of the merons and antimerons. In this respect 
it is worth noting that dissipationless spin current generation was recently predicted and 
observed in systems with spin-orbit interaction , 9 i 10 i n where the spin current is produced by 
a spin Berry phased induced by spin-orbit interaction. 

For the appearance of a persistent current, however, one more twist is necessary. Al- 
though merons and antimerons generate a spin current, electric current is absent due to 
the exact cancellation of the up- and down-spin electron currents. This is an example of 
the so-called Block's theorem which may be stated: " The ground state of a time- reversal 
invariant Hamiltonian is currentless".— In the present work we argue that the Bloch's theo- 
rem is circumvented through the fact that states with merons and antimerons are stabilized 
by an imaginary magnetic flux that "negates" the fictitious magnetic flux from the Berry 
phase; namely, if a flux that counteracts the distructive interference of the wave function 
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around each core of meron or antimeron caused by the spin Berry phase (note that the Berry 
phase causes destructive interference due to 7r flux through the cores), the total energy is 
significantly decreased. Although the introduction of such an imaginary flux is artificial, it 
becomes a real one when the small polarons form side-by-side bipolarons; this realization of 
the imaginary flux occurs because the vector potential for the imaginary flux which is writ- 
ten as Aflux = 7px( r ) (<? is electron charge and x is an angle- variable with periodicity 2ir), 
can be transfered from the Hamiltonian to a phase factor on the wave function written as 
exp(— i| x( r j))i where N is the total number of electrons and is the coordinate of the 
jth electron. In other words, the Bloch's theorem is circumvented by the spontaneous gen- 
eration of the motion of the system as a whole given by a wave function exp(— i| £^ x( r j)) 
by itself. Note that this current carrying state has a partner state that is connected by the 
time-reversal operation; thus, the stability of the persistent current state depends on the 
tunneling life time to the partner state through the energy barrier between them. 

An interesting effect is predicted for the above current carrying state if the bipolarons 
move; such a collective mode will give rise to a time-dependence of x- Since time variation 
of a phase of a wave function gives rise to a chemical potential shift,— which for charged 
particles means an appearance of an electrostatic potential whose electric field is given by 
E = — ^ dA g t u * = — I^Vx- Then, the Nernst effect is expected by the thermal diffusion of 
the bipolaronic cores of merons and antimerons. This may explain the anomalously large 
Nernst signals observed in the pseudogap phase of the cuprates.— 

As will be shown later, the present model also gives qualitative explanations for the 
appearance of the Drude-like peak in the optical conductivity measurement and the Fermi 
arc observed in the angle-resolved photoemission spectra in the cuprates; and suggests the 
appearance of the phase variable or "Nambu-Goldstone" mode for the superconductivity 
without the breakdown of the global U(l) gauge invariance in the cuprates. 

This paper is organized as follows: in the next section we discuss merons and antimerons 
embedded in a spin density wave background in a two dimensional square lattice using a 
Hubbard model. In Section III, it will be shown that a Berry phase arises from the spin 
rotation of conduction electrons in the spiral spin texture of merons and antimerons, and 
generates a spin current. In Section IV, we will add an imaginary magnetic flux that negates 
the fictitious magnetic flux from the spin Berry phase, and show that such a flux significantly 
lowers the total energy of the state with merons and antimerons. In Section V, we point out 
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that imaginary flux added in Section IV becomes a real one when side-by-side bipolarons are 
formed. In this case, the flux arises spontaneously from a wave function for the system as a 
whole, and circular electric currents around the bipolarons are produced. The appearance 
of the Drude weight is also expected in such a situation. In Section VI, we discuss the 
appearance of large Nernst signals by the flow of bipolarons. Lastly, in Section VII, the 
relevance and implications of the present results in cuprate superconductivity are discussed. 
The present work suggests a "Nambu-Goldstone" mode appears without the breakdown of 
the global U(l) gauge invariance in the cuprates. 

II. MERONS AND ANTIMERONS IN A SPIN-DENSITY WAVE BACKGROUND 

The model we are going to use is a Hubbard model with including small polaron effects 
of doped holes. The Hubbard Hamiltonian we use is given by 

where c\ a (c ia ) denotes a creation (annihilation) operator for electron with spin a at site i, 
and the first sum is taken over pairs of sites. 

We employ a mean field approximation that gives rise to a spin density wave (SDW) 
ordered electronic state with spin polarization in the x-y plane. Such a spin modulation 
is relevant for 2D organic conductors,— and also to the cuprates.— Electronic states are 
obtained from the following mean field Hamiltonian derived from Eq. @: 

H **f = - + h -c-) - u ( A /:)- r /. + A /:i. r r) 

{i,j),cr 3 

+ u \ a j\ 2 + 1 E(sV;t + c] lCjl ) (6) 

3 j 

where the mean field at site j is given by 

A, = |A,|e^ = (c],c n ). (7) 

When the Coulomb interaction parameter U is much larger than any of the absolute 
values of transfer integrals \tji\, the second sum in Eq. (JBJ) is the dominant term. In this 
situation, the representation that diagonalizes the second term is useful, 

- U E ( A i c ,W + = -u E lA^I - > ( 8 ) 

3 j 
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where the new annihilation operators aj and bj are related to the original by 
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Then, the first sum in Eq. © becomes 



- ^ k cos ~ 2 ^ ( a 3 afc + h ) hk > + il i k sin ^ 3 - — (a)b k + b]a k ) + h.c.]. (10) 
At half filling in the large U limit, the ground state is given by 

i°> = n«ji vac >- ( n ) 

3 

Now, we take into account the first nearest neighbor hopping term, 

-it ^2 sin ^%4 6 * + &K) + h - c -> ( 12 ) 

where a pair of first nearest neighbor sites k and j is denoted as (k,j)i. Using the pertur- 
bation theory, we obtain the following minimal energy condition for the phase £: 

sin 2 = i. (13) 

This condition is satisfied when the phase £ at site j is chosen as 

C, ~(.i.r ■ jy), (14) 

where (j x ,j y ) is the coordinate of the site j, and the lattice constant is set to unity. The 
above phase yields an antiferromagnetically ordered state in accord with the well-known 
result at half filling. Note that the antiferromagnetic order here are not due to the Fermi 
surface nesting; it is due to the strong Coulomb interaction augmented by the first nearest 
neighbor hopping perturbation. Therefore, the mean field treatment employed here captures 
the physics of the Mott insulator to some extent. 

Taking into account the antiferromagnetic correlation, the energy dispersion for the lower 
Hubbard band is calculated as 

e(k) = ^ + Q) + e(k) _i^ (k + Q) _ f(k)|2 + t/2| (l5) 

where Q = (tt, it) and the single particle energy e(k) is given by 

e(k) = — 2t(cos k x + cos k y ) — At' cos k x cos k y + 0.5U, (16) 
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FIG. 1: Energy dispersion of single particle states for an antiferromagnetic spin density wave state 
e(k). The parameter used are t'/t = —0.2 and U/t = 8 with t = 1. The lattice constant is taken 
to be unity 

with t and t' being the first and second nearest neighbor transfer integrals, respectively (the 
first and second nearest hopping terms are included in this work). The parameters t, t', 
and U are chosen to be t'/t = —0.2 and U/t = 8. With these parameters the dispersion 
in Eq. (1 15 J) resembles the one seen in a parent compound of a cuprate superconductor. In 
Fig. [IJ the single particle energy dispersion e(k) is depicted. It shows a similar dispersion 
observed in the ARPES measurement for Ca2Cu02Cl2. 1 ^ 

When holes are doped, we assume that they become small polarons due to strong electron- 
phonon interaction and the strong Coulomb interaction, which seems to be relevant to the 
cuprates . 2 i 3 i 4 i 5 i 6 i 19 Further, we add the condition that the hopping motion of the polarons is 
so slow that polarons can be treated as static objects. In this situation, a part of the small 
polaron effects can be taken into account simply by setting the transfer integrals from and 
to polaron sites to zero. We will use this static approximation in the following. Furthermore, 
we assume that the site energy for the polaron occupied site is such that the small polaron 
formation is energetically favorable; i.e., the polaron level is above the upper edge of the 
filled lower band. We achieve this by adding the term E s (c^c^ + c^Cji) in Eq. ©, where 
E s is the polaron stabilization energy stems from the vibronic interaction.— For simplicity, 
we take E s = — \U in the following calculations; namely, we use the mean field Hamiltonian 
given by 

H 'mf = ~ E M<, c - + h.c) -UJ2 {^W + \'i'\r) +UJ2\^\ 2 - (17) 

(ij),o- j 3 

The electronic state is calculated using the mean-field Hamiltonian Eq. ()17j) by self- 
consistently obtaining Aj using Eq. ((7j). The number of electrons is the number of sites 
minus the number of polarons (or holes). In Fig. Efa), the spin texture for the state with 
the phase £j = 7i(j x + j y ) is depicted; the system size is a 16 x 16 lattice, and the periodic 
boundary condition is employed. 

Now, we consider the case where doped holes become cores of merons and antimerons. 
In general, merons and antimerons appear in meron-antimeron pairs because such paris do 
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FIG. 2: Spin textures for states with four small polarons embedded in an antiferromagnetic 
background. The four small polarons are located at (6, 6), (9, 6), (6, 9), and (9, 9). (a): spin texture 
for the state without a meron or an antimeron. The total energy is — 614. 5t; (b): the same as (a) 
but the polarons are now cores of merons (indicated by M), and antimerons (indicated by A). The 
total energy is higher than that for (a) by 0.7t. 

not severely destroy the background spin texture. Thus, the conservation of the winding 
number is expected since the winding number for a meron is +1 and that for an antimeron 
is -l.fi 

In order to include the spin texture modulation by merons and antimerons, we write the 
phase £ as 

& = *u*+Jv) + E tan " - E tan_1 ( 18 ) 

^ Jx-Xm ^ 3x-x A 

where M and A in the sums run through merons and antimerons, respectively. In the present 
work the phase £j is fixed; we tried the optimization for £j by the procedure given in RefJ^, 
but it did not lower the total energy in the present mean field treatment. 

In Fig. E(b), the spin texture for the state with two merons and two antimerons is shown. 
The total energy for this state is — 613. 8£, while the one without a meron or an antimeron is 
— 614.51 This shows that the formation of merons and antimerons destabilizes the system. 
This contradicts the results in Refl^f] where the stabilization by the meron formation is 
observed; this difference probably comes from the neglect of the nearest neighbor Coulomb 
interaction and the inclusion of the small polaron effects in the present calculation. 

III. SPIN CURRENTS GENERATION BY MERONS AND ANTIMERONS 

In this section, we show that circular spin currents appear around polaronic cores of 
merons and antimerons. 
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FIG. 3: Currents for states with the spin texture depicted in Fig. |2fb). Three different cases with 
and without an imaginary flux through cores of merons and antimerons are shown, (a): current 
for up-spin electrons for the case without an imaginary flux. In this case the current for down-spin 
electrons is exactly opposite to that for up-spin electrons. The total energy is — 613. 8i; (b): the 
total current for the state with Afj ux i. The total energy is lower than that for (a) by 2.6i; (c): 
the same as (b) but with Afl UX 2. The total energy is lower than that for (a) by 3.0t. The electric 
currents are obtained by the multiplication of — e. 

Current densities are calculated using the following current density operators; 

s 

+ {^l+x+ys C ls ~ c \ s c l+x+ys + c\+x~ys C ls ~ c \ s C l+x-ys^ > 

* s 

j y (l) =itJ2 ( c \ +y s c is - 4s c i+ys) 

s 

+ Vf*' ^ ( C WW* - C lci + x + ys + 4_ x+ys C ls - Ctc^ x+ys ) , (19) 

where I + x, I + y, I — x, and I — y denote the nearest neighbor sites of I in the x, y, —x, 
and — y directions, respectively; and I + x + y, I — x + y, I + x — y, and I — x — y denote 
the second nearest neighbor sites of I in the x + y, —x + y, x — y, and —x — y directions, 
respectively. Note that the above operators may be deduced from the Heisenberg equation 
for cj^Cjcr.— The electric current densities are obtained by the multiplication of the charge 
q = — e on the current densities, where — e is the electron charge. 

The electric current densities calculated for the states with spin textures depicted in 
Figs. Ufa) and (b) are zero. This is in accord with the Bloch's theorem, although the latter 
has nonzero spin current. In Fig. 3 (a), the current for up-spin electrons is depicted (the 
current for down-spin electrons is exactly opposite). 

The origin of the spin current can be attributed to the phase factors e l ^ J and e~ % ^ j 
in Eq. (J5jl. which become multi- valued when merons and antimerons are present. A multi- 
valued phase factor (or a non-integrable phase factor) is a hallmark of a "magnetic field" ,— 
and such a multi-valuedness can be ascribed to a fictitious magnetic field. Let us find 
the vector potential for the fictitious magnetic field. We denote Wannier functions for the 



annihilation operator cj s by Wj. We see that when the new operators, aj and bj in Eq. Q, 
are employed, the phase factors e l ^ j and e~ l ^ j are multiplied on and Cjj, respectively; 
this is interpreted that Wannier functions are modified as 

w n = Wj e^, (20) 

for Cji, and 

w n = Wje-ifc (21) 

for Cji- According to the prescription for the Peierls substitution , 23 ^ 4 the above modifications 
of the Wannier functions are equivalent to those performed when a system is in a magnetic 
field with vector potentials 

ch 

A/ lcT (x i ) = -^V£„ (22) 



for up-spin electrons, and 



A/fcife) = -V0, (23) 



for down-spin electrons. The current for up-spin electrons is due to the fictitious vector 
potential Ay^-p, and that for down-spin electrons to Af ic f, they are exactly opposite, Aj ic | = 
—Af ic i, and produce opposite currents, jj = — 

Strictly speaking, the transformation matrix in Eq. Q is not single- valued when merons 
and antimerons exist; thus, it is not a legitimate unitary transformation. Actually, if we 
follow the standard calculational procedure for the Berry phase, we obtain the same fictitious 
vector potentials with using a legitimate unitary transforamtion as follows: 25 ' 2 ^ we use the 
following single-valued transformation matrix instead of the one in Eq. Q 

V * = -ft \ , '( ii£ ) ' (24) 

where n and m are integers, and Vj is a single-valued unitary matrix. 

Associated with the above unitary matrix is a 2 x 2 Berry connection given by 

-iV j VV=-^\ . (25) 

2 I ne <m-n) 2m -1 
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In the following, we employ the adiabatic approximation where off-diagonal elements of 
the above Berry connection that causes transitions between the lower and upper bands are 
neglected. In the situation where the energy gap between the upper and lower Hubbard 
bands is much larger than external perturbations and the upper Hubbard band is empty, 
this approximation is a good one. Using this approximation, the motions of electrons in the 
lower Hubbard band are separately treated from those in the upper one. 

From Eq. (|24|). the annihilation operator for electrons in the lower band at site j is given 

by 

a^-L^c^ + e^-^c^), (26) 

and from the (2, 2)-element of the Berry connection in Eq. (|25|) the fictitious vector potential 
for this branch is obtained as 

(2m-l)^V£ r (27) 

For up-spin electrons, an additional vector potential -m^V(j arises from the phase 
factor e m ^ J on Cjf in Eq. (f2T)|) . The total vector potential for the up-spin electrons is the 
sum of Eq. (j2T|) and — m— V£j, i.e., — §^V£j, which agrees with Eq. (j22J)- Likewise, for 
down-spin electrons, an additional vector potential — (m — l)*V£j arises from the phase 
factor e^" 1-1 ^ on Cji in a'j. The total fictitious vector potential for down-spin electrons is 
given as the sum of Eq. (|2*7)l and — (m — 1)— V£j, i.e., j*V£ 3 - in agreement with Eq. (|2*3Jl . 



IV. IMAGINARY FLUX STATES 



In the previous section, it is shown that a spin current is produced by merons and 
antimerons; but the electric current is absent because the up- and down-spin currents exactly 
cancel each other. In this section, we add an "imaginary magnetic flux through cores 
of merons and antimerons. If we add a suitable flux it is show that states with merons 
and antimerons are stabilized, and gives rise to non-zero electric current; however, we do 
not concern how such a flux arises in reality in this section (for this reason we call it an 
"imaginartf flux); the mechanism for the appearance of such a flux will be discussed in the 
next section. 
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Let us apply a magnetic flux given by the vector potential 



A fluxl (x) = ^Vxi(x), (28) 



where xi is given by 



Xl (x) = £tan-'^M-l>n->^i. (29) 
x-x M ~ x-x A 

The flux Eq. (|28|) is chosen since it may "negate" the fictitious magnetic field produced by 
merons and antimerons, and stabilizes the states with merons and antimerons. The negation 
here does not mean total cancelation of the effect of the fictitious magnetic field, but means 
the retrieval of the constructive interference of the wave functions around the cores from 
the destructive interference caused by the fictitious magnetic field from the Berry phase. 
The negation is expected because 1) the vector potential A fluxl cancels the fictitious vector 
potential for the up-spin electrons, 2) the sum of A fluxl and Af ic ^ only gives rise to an 
integrable U(l) phase factor for the down-spin electrons. 

In Fig. 0(b), the electric current for the state with A fluxl is depicted. Currents are 
calculated using Eq. (JT5J) but modifying the transfer integrals according to the following 
Peierls substitution^^ 

f - f pxn (a A fluxl ( rj ) + A fluxl (r fc ) \ 

tjk - t jk ex P \ n ^ l r j - r k) J ■ (ouj 

Note that we will use a similar modification of the transfer integrals when we calculate 
current densities for states with a flux different from Afl ux i. Significant stabilization is 
brought about by the introduction of the imaginary flux; the total energy is reduced by 2.6t 
for the case seen in Fig. 0(b), which is even lower than that for the state without a merons 
or an antimeron by 1.9t. 

Actually, many other choices are possible for the imaginary flux that negates the fictitious 
magnetic field. Next we try the following vector potential 

cfi 

A flux2 (x) = — Vx 2 (x), (31) 
Zq 

where \2 is given by 

Wx) = } tan h y tan . (32) 

/ — ' X — Xm X — Xa 

M A A 
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In Fig. Efc), the current for the state with Aa UX 2 is depicted. The current for this state is 
much larger than the one with Afl ux i. Its total energy is — 616. 8t, which is lower than that 
for the state with A fluxl by 0.4t. 

In order to assess the stability brought about by the introduction of the imaginary flux, 
energies for three different cases with two polarons in the 16 x 16 lattice are calculated: i) 
case without a meron or an antimeron; ii) case with a meron and an antimeron but without 
imaginary flux; iii) case with a meron, an antimeron, and an imaginary flux. The results are 
listed in Table |l] For the most of the cases, the state with a meron, an antimeron, and the 
vector potential Afl UX 2 is the lowest in energy among the cases considered. Note that the 
energy for states without a meron or an antimeron is almost always — 620. 6t (or 0.3t above 
the lowest one in Table H}. The state becomes more stable when the meron and antimeron sit 
closer in accord with the well-known fact that attractive force exists between a meron and an 
antimeron.^ The most stable state is the one with an antimeron at (7, 6), i.e., the state with 
the side- by-side meron-antimeron pair, and Afl UX 2. As will be discussed in the next section, 
a mechanism exists that gives rise to the imaginary flux in this side-by-side bipolaron case; 
thus, this indicates the spontaneous formation of side-by-side bipolarons with Afl UX 2. It is 
also worth noting that between the two types of the vector potentials Afl uxl and Afl ux2 , the 
latter always gives slightly lower energy. 



V. CURRENT CARRYING STABLE STATES 



In this section, we show that the imaginary flux added in the previous section becomes a 
real one when side-by-side bipolarons are formed. 

Let us take a side-by-side bipolaronic core of a meron-antimeron pair, and consider the 
gauge potential for the imaginary flux for it, 



A p flux2 = !|V;& (33) 



where the angle variable x p is given by 



X§(x) = tan- 1 + tan- 1 ^— ^; (34) 



ry ry 

' M A 



ry ry 1 ry ry- 1 

i x Mi Vm) an d i x A.y Va) denote locations of the cores of the meron and antimeron, respectively. 
An example of electric current produced in the state with a side-by-side bipolaron and Ag ux2 
is depicted in Fig. HJ 
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TABLE I: Total energies (in the units of t) obtained by the mean-field Hamiltonian. The energy is 
relative to the lowest one, i.e., the case for (6, 7) with a meron, an antimeron, and Afl ux2 . In each 
box, total energies for a case without a meron (M) or an antimeron (A), and that with a meron 
and an antimeron (M and A with Afl ux2 ) are listed. One of the two polarons is always placed 
at (6,6) and the other at (j x ,jy) in the 16 x 16 square lattice. The meron core is always placed 
at (6,6) and the antimeron core at (j x ,jy)- For the side-by-side meron-antimeron core cases, the 
energy with Afl uxl is also given. 
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3.0 
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M and A (A fluxl ; A flux2 ) 


3.9(1.2) 


4.2(1.6) 


4.5(2.0) 


4.9(2.3) 


jy = 7 


3.0 


3.0 


3.0 


3.0 


3.0 




3.3 (0.9; 0.3) 


3.7(1.0) 


4.1(1.3) 


4.4(1.8) 


4.8(2.1) 


jy = 6 


2.8 


3.0 


3.0 


3.0 


3.0 




2.9 (0.2; 0) 


3.6(0.8) 


4.0(1.1) 


4.4(1.6) 


4.7(2.0) 



For every loop delineated by the hopping motion of electrons, the change of |x 2 * s a 
multiple of 2ir since it encircles both of the cores, or neither of them (note that the transfer 

■ 1 V N p ( \ 

integrals from and to polaron sites are set to zero); then, the phase factor e - ^ 2 ^' x ^> 
becomes single-valued as far as the conduction electron motion is concerend. Therefore it is 
legitimate to transfer the vector potential from the Hamiltonian to a wave function as 

tt'Cxx,... , XJV ) = e -^fx^')*( Xl ,... , XJV ), (35) 

where the is the resulting single-valued wave function and \l/ denotes the original one. 
Actually, the phase factor e~* 52 ^ x 2 i x jJ j s re g ar ded as a wave function for the motion of the 
system as a whole. This whole system motion will arise spontaneously since the state with 
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FIG. 4: Current for the state with side- by-side cores of a meron and an antimeron, and Afl ux 2. 
The meron core is located at (6, 6) and the antimeron core is located at (7, 6) of the 16 x 16 square 
lattice in Fig. [2J The electric current is obtained by the multiplication of — e. 



FIG. 5: Currents generated by side-by-side meron-antimeron cores, and the flux Ag ux . (a): 
(a&,!&)=(6,6), (10,6), (7,10), (12,10) and {x%,y%)=(7,6), (11,6), (6,10), (11,10). The total 
energy is E = -606.58*; (b): (a&,!/&)=(6,6), (7,10), (a£,l/£)=(7,6), (6,10), (x%, |&)=(10, 6), 
(12,10), and (x% I/J)=(ll, 6), (11,10). The total energy is E = -606.08i; (c): (x^, y^)=(4, 2), 
(5,6), (4,10), (5,14), (x^,y^)=(5,2), (4,6), (5,10), (4,14), (a&, y£)=(12, 2), (13,6), (12,10), 
(13,14), and (sj, yj)=(13, 2), (12,6), (13,10), (12,14). The electric currents are obtained by the 
multiplication of — e. 

it is lower in energy than that without it. 

Let us consider the case where many side-by-side meron-antimeron cores exist; we adopt 
the gauge potential for the flux written as 

cfi 

A flux (x) = ^Vx(x), (36) 

where the phase x is given by 
X (x) = y( tan- 1 + tan" 1 - V [ tan' 1 ^| + tan" 1 ) ;(37) 

the first sum corresponds to contributions from pairs with the gauge potential analogous 
to 

Afl UX 2 that produce clockwise currents around merons and antimerons; and the second 
sum to — Ag ux2 that produce counterclockwise currents around merons and antimerons. In 
Fig. El examples of electric currents generated by various Afj ux are depicted. As seen in 
Fig. an extended current is generated by a collection of side-by-side bipolarons. The 
total energy for the state in Fig. E[a) is slightly lower than that for Fig. Efb); this seems to 
indicate that a larger circular current state is energetically favorable if the magnetic energy 
is neglected. If the magnetic energy is included a large single domain circular current state is 
not energetically favorable; current domains are formed to minimize the sum of the magnetic 
energy and total electronic energy. 
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FIG. 6: Single particle energies for the state with the current shown in Fig. Etc). The unit of energy 
is t. The levels in the lower Hubbard band are fully occupied, and those in the upper Hubbard 
band and in-gap states are empty. The in-gap states are due to small polarons whose position 
depends on the polaron stabilization energy E s ; in the present work it is chosen as E s = —\U. 

Note that the current obtained here is the one outside of the conventional band theory; 
the charge carriers are electrons in the fully occupied lower band. Fig. |U] shows the plot of 
single particle energies for the state with the current shown in Fig. |3^c). The lower Hubbard 
band is fully occupied, yet, produces the electric current. 

The Drude weight is also expected in this system: the total energy is a functional of 
Aflux, written as E tot [ Afl ux ] ; the total energy depends on Afl ux , and should be bounded 
below, thus, an optimal vector potential A^ ux should exist; then, a general vector potential 
A flux is written as a sum of the optimized one and a deviation from it, A flux = Ag ux + 5Afl ux . 
Assuming that the deviation 5Afl ux is small compared with the optimized one, the total 
energy is written up to the second order in 5Afl ux as 

£ tot [A flux ] w E tot [A° Qm ] + / d 3 x d 3 y *(5A flux (x)) 5A Rux (y). (38) 

J dA flux (x)dA flux (y) 

Thus, the Drude weight (in the x direction) is calculated as^ 



7TC 



2' 



d E tot [A Rux + A x e x ) _ 2 f 3 3 S EtotfA^ 



- vrc 2 / d 3 x d'y — \ ;TJT\ < ^ > ( 39 ) 
A x =o J o(A flux ) :E (x)d(v4 flux ) :r (y) 



dA 2 x 

where is the unit vector in the x direction. The Drude weight describes the sensitivity of 
the system with twisting the periodic boundary condition; if a collection of circular currents 
around the bipolarons extends over the sample, the nonzero Drude weight is expected. 

In the current carrying state considered in this section, the Bloch's theorem is circum- 
vented by the spontaneous generation of the whole system motion by itself; it produces a 
vector potential that negates the fictitious magnetic field from the spin Berry phase, and 
stabilizes the state with bipolaronic cores of merons and antimerons. Since the Hamiltonian 
is time-reversal invaraint, this current carrying state has a partner state that is connected 
by the time-reversal operation. These two states may be viewed as states at two minima 
of a double-well potential. A current carrying state is stable if the energy barrier between 
them is so large that the tunneling motion between them is practically negligible. 

16 



Lastly in this section, we briefly discuss the response of the system to an external magnetic 
field. The state with Afl ux and that with — Afl ux have the same energy and produce mutually 
opposite currents; when an magnetic field is applied in the z direction, their energies become 
different. Let us apply a magnetic flux density Be z in the z direction to the state (e 2 is the 
unit vector in the z direction) with the current shown in Fig. |Sfa); the result is the total 
energy becomes — 606. 55t for B = —0.001, and — 606. 47t for B = 0.001 where B is given in 
the units of 2$o/a 2 with a being the lattice constant and 2$o equal to he/ e=2.07 x 10~' 
gauss-cm 2 . This result shows that the one with the diamagnetic current to the external 
magnetic field is more stable (the energy difference 0.08t is about 200 K for the cuprates). 
In many cases (but not always) we have examined, the state with a diamagnetic current 
seems to become slightly lower than that with a paramagnetic current; this may imply an 
appearance of a large diamagnetic response. 



VI. NERNST SIGNALS BY THE FLOW OF MERONS AND ANTIMERONS 

When the flow of bipolaronic cores of merons and antimerons occurs, the phase x becomes 
time-dependent. This means that conduction electrons are under the influence of an electric 
field given by 

l<9A flux h 

E = --— = --V* (40) 

In this section we show that this electric field gives rise to Nernst signals. 

To calculate the electric field produced by a thermal gradient — VT, let us take a rectan- 
gular system shown in Fig. [7|with the length L x in the x direction (0 < x < L x ), and the 
width L y in the y direction (0 < y < L y )\ a magnetic field is applied in the z direction; and 
the thermal gradient is in the x direction —d x T. 

From Eq. P0j) . the voltage across the width L y measured at x = is calculated as 

h f Ly , d (L T 



V y = E y L y = — I dy—x [ —,y 



2eJ »dy^\2 



27i 



(41) 



After the movement of a bipolaron from x = to x = L x , the difference of the phase x 
across the width at x = % changes. If L x is sufficiently large, the phase change is — 2n if x 
is the one produces a clockwise circular current (27r for the one produces a counterclockwise 
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FIG. 7: Experimental set-up for the Nernst signal measurement. 



circular current) around a meron or an antimeron. In the situation where the bipolaron 
dissociation and recombination processes are frequent, the above phase change should be 
regarded as an average over many disconnected flows of bipolarons. Similarly, in an average 
sense, we denote the velocity of the flow by v, and use L x /v as the average time duration 
for the movement of a bipolaron from x = to x = L x . 

Then, the voltage measured in the Nernst experiment is calculated as 

V - (42) 
y ~ 2e L x /v ' 1 ' 

where N m is the difference of the numbers for diamagnetic circular current cores (No) and 
paramagnetic circular current cores (N p ) of polarons, N m = N d — N p . N m is positive in the 
presence of the applied magnetic field because diamagnetic current states are expected to 
be more stable than paramagnetic current states. 
Finally, the electric field produced is calculated as 

E _ Vy_ _ hv N m = hvn m ^ 

where n m is the surface density of N m given by N m /L x L y . Since a loop current gives 
rise to a magnetic moment, circular currents around the bipolarons produce magnetization 
perpendicular to the two dimensional plane. If we assume that each meron or antimeron 
core in a bipolaron contributes to roughly the same magnetic moment, Eq. (J43|) indicates 
that the Nernst signal is roughly proportional to the magnetization. Actually, the Nernst 
signal proportional to the magnetization has been observed in the pseudogap phase of the 
cuprate superconductors.— 



VII. DISCUSSION 



In this section, we shall discuss the relevance and implications of the present work in the 
cuprate superconductivity. 

We have been dealing with the inhomogeneous system composed of conduction electrons 
and small polarons; the former constitute a background SDW order, and the latter become 
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embedded bipolaronic cores of merons and antimerons. For lightly doped cuprates, the 
strong evidence of polaron formation exists ; 28 > 29 and the evidence of polaronic charge carriers 
also exits4^ 

The single particle levels depicted in Fig. |H1 show the splitting of the original band into 
the lower and upper Hubbard bands, and between them are the in-gap states due to po- 
larons. In the lightly doped cuprates, the constant chemical potential with changing the 
hole concentration has been observed; it is interpreted due to the presence of in-gap states 
that pin the Fermi level . 31 ^ 32 Spectral weight transfer that appears to be between the in-gap 
states and the lower Hubbard band has been observed.— Quite often, the charged stripes 
are argued to be responsible for the chemical potential pinning;— and the splitting of peaks 
in the neutron scattering experiments has been taken as the evidence for the existence of 
them . 35 i 36 In contrast to it the present work suggests that the chemical potential pinning 
is atributed to the small polaron formation; those polarons form merons and antimerons 
around them and explain the neutron scattering peak splitting; 8 thus, what has been con- 
sidered to be due to charged stripes may be actually due to merons and antimerons whose 
cores are small polarons. 

It has been though that a disconnected Fermi surface like feature (called "Fermi arc") 39l4 ° 
seen in the angle resolved photoemission spectroscopy (ARPES) and the Drude-like peak 
observed in the optical conductivity measuremen t 37 ! 38 indicate that the charge carriers are 
not small polarons, thus, small polarons are thought to be not relevant to cupreate super- 
conductivity. However, the present model explains those phenomena with the small polaron 
as an essential ingredient, although it is not the charge carrier. 

Let us first simulate the ARPES experiment: we calculate the ARPES plot as a sum 
of the square of the absolute value of the Fourier transform of single particle level wave 
functions given by 



where %\ is the arbitrarily chosen lower limit of the occupied state, i% is the highest occupied 
state in the lower Hubbard band, and (pi is the ith single particle wave function. The results 
are depicted in Fig. |H1 A Fermi arc which is similar to the one observed in ARPES^ is 
seen in Fig. |H] The Fermi arc is more arc-like when the circular currents are all in the same 
(clockwise or counterclockwise) direction and produce a large current region. Therefore, the 




(44) 
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FIG. 8: Simulations of the angle-resolved photoemission spectra (ARPES) for cuprates. Each row 
has three plots; Spin texture plot, the total current plot, and ARPES simulation plot (symmetrized 
about the (0,0) — (vr, 7r) reflection line), from left to right. The first and second rows are those for 
states with one bipolaronic core. The ARPES plot is calculated using Eq. (|44|) with the 15 highest 
occupied single particle level wave functions (energy range is from — 3.92t to — 3.80i for the first 
row and — 3.93i to — 3.80t for the second). The total energy for the first row case is lower than that 
for the second row case by 0.2t. The third row is for the state with 16 bipolaronic cores. In this 
case, the Coulomb interaction parameter is reduced from U = 8t to U = 4t; this change causes a 
slight elongation of the arc compared with the U = 8t case. The ARPES plot is obtained using 
Eq. ()44|) with the 25 highest occupied single particle level wave functions (energy range is from 
-1.90* to -lMt). 

experimentally obtained Fermi arc suggests the presence of bipolaronic cores of merons and 
antimerons with a large circular current area, which is in accord with the result in Section 
IV that a larger current state is energetically favorable. 

Next, we shall consider the Drude peak problem. The Drude weight of Eq. is expected 
to arise in states with circular currents around small polarons. Low energy excitations 
associated with nipping of circular currents (from clockwise one to counterclockwise one, or 
vise versa) exist; an estimate from the case depicted in Fig.lHlis that it is of order of 0.2t « 40 
meV. Therefore, the Drude-like peak observed in the optical conductivity measuremen t'^ 1 ^ 
is qualitatively explained in the present model. 

In the cuprates, the anomalously large Nernst signal and enhanced diamagnetism have 
been observed at temperatures well above T c . 15 i 41 A large Nernst signal is usually associated 
with the flow of Abrikosov vortices, thus, it is tempting to regard the enhanced Nernst 
signal as evidence for the spontaneous formation of moving Abrikosov vortices. Besides, 
it is also in accord with the scenario that the onset of superconductivity corresponds to 
the establishment of a long-range coherence of the order parameter phase. 42 However, it has 
been argued that the presence of a supercurrent that is capable of forming Abrikosov vortices 
seems to be implausible at such high temperatures since the sharp resistive and magnetic 
transitions at T c cannot be explained by the fluctuation scenario.— Indeed, the scanning 
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FIG. 9: Temperature dependence of the magnetization for the x = 0.16 (Bi 2212) sample obtained 
by Y. Wang et al»^ Filled circles are experimental data. The fitting function is M(T) = —Ae~ T ^ T 
with A = 0.00033 A/m and T f = -1, 100 K. 

SQUID microscopy observed the disappearance of the Abrikosov vortices above T c .— The 
origin of the large Nernst effect is a very controversial issued The present work provides a 
new explanation: it is due to the flow of the bipolaronic cores of merons and antimerons. 
The comparison of the experimental results 15 and the formula in Eq. indicates that the 
magnitude of the magnetization increases exponentially with the increase of 1/T, which may 
be attributed to the increase of the number of the bipolarons; if we assume each bipolaron 
contributes either an average up or down magnetic moment produced by the circular current 
around it, the magnetization will be proportional to the difference of the numbers of the up 
and down magnetic moments per volume. Let us try a fit to the experimentally measured 
magnetization to the function M{T) = — Ae~ Tf / T , where A and Tf are fitting parameters; 
the fit is very good as shown in Fig. 121 From the fit, the formation energy —Tf of the 
bipolaron is estimated to be 75 meV, which is close to the value where a kink is observed 
in the energy dispersion of angle resolved photoemission spectrum.^ Note that this energy 
seems to be larger than that is expected for the flipping of the circular currents; and may 
correspond to the formation energy of a bipolaron from two separate polarons. 

In terms of the present model, the experimental fact that the observed Nernst signal above 
T c smoothly connects to that by the Abrikosov vortices below T<J^ implies a rather surprising 
conclusion: the flow of the bipolarons smoothly connects to the Abrikosov vortex flow, thus, 
the phase that generates the supercurrent below T c (the phase of the order parameter) is 
actually the phase that generates a circular current around the bipolaron above T c (the 
phase for the gauge potential A flux ). In other words, the cuprate superconductivity may be 
the superconductivity with the bipolarons serving not as the charge carriers, but as pumps 
for the supercurrent. 

One may think that in such a theory, the flux quantization in the units of or the ac 
Josephson effect with the frequency 2eV/h (V is the applied voltage) will not be possible 
because the charge on the current carrier is not 2e. However, as will be explained below, the 
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present theory will explain the observed flux quantization and Josephson frequency. Note 
that 2e in ^ or 2eV/ h does not exclusively mean the existence of 2e charge carriers.— 

Let us examine the flux quantization problem: we first show that when an external 
magnetic field given by the vector potential A em is added, the sum A em -|- Afj ux is gauge 
invariant. For the gauge transformation 

A cm -> A cm + V0, (45) 
the total wave function transforms as 

e -*3E i= ix(ri)^( ri) ... , rff ) e »E J =ine L e -»iE,= 1 x('i)<f ( rij . . . >rjv ) j ( 46 ) 

thus, the phase \ transforms as 

_ 2# 

X ¥ X »- > 
he 

which leads to the transformation of Afl ux as 

A flux -> A flux - V0. (48) 

Combining Eqs. ()45|) and (HHJ), it is deduced that the sum A em + A flux is gauge invariant; in 
other words, the phase \x is the superpotential of F. London 48 that absorbs the arbitrariness 
in the gauge. 

We can calculate the electric current using the energy functional in Eq. (|38|) by substi- 
tuting A flux + A em for A flux ; 

£tot[A flux + A em ] sa E tot [A° aux } 

+ \ [ d 3 x d 3 y *(5A flux (x) + A cm (x)) -J^^d. (5A flux (y) + A cm (y)) , (49) 
J c'A flux (x)dA flux (y) 

where <5Afl ux + A em is regarded as a small perturbation. The expansion above usually does 
not exist if the gauge potential A flux is absent because the arbitrariness in the gauge of A em 
makes such an expansion meaningless; however, thanks to the presence of Afl ux , the total 
energy is now a functional of the gauge invariant sum Afl ux + A om , thus, the above expansion 
does not depend on the choice of the gauge and is meaningful. In Eq. (J3HJ), the term 5Afl ux 
is added to take into account the slight change of A flux from Ag ux due to the presence of the 
magnetic field V x A em ; thereby, the expansion becomes the one in terms of S Afl ux + A em . 
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(47) 



Using the functional Eq. (fl§|). the electric current is calculated as 

J(x) = -c[ d 3 y ^f E ;°\ [ ^ J ( ; (5A flux (y) + A cm (y)) . (50) 

This is gauge invariant since 5Afl ux + A cm is gauge invariant. It is known that the current 
given above leads to the Meissner effect, which is in accord with the fact mentioned in 
Section V that the state with the diamagnetic current becomes more stable when an external 
magnetic field is applied. 

The Meissner effect means that the current is zero deep inside a large enough sample. 
Then, we have the relation, A om = — <5A flux = — which leads to the flux quantization 

$ = j)A cm -dr = -^(j)5\7x-dr = 0ra, (51) 

where n is an integer. In order to have a nonzero n, the conservation of the winding number 
must be lifted. Actually, when the quantized flux |^ occurs, the spin texture is the one 
produced by an isolated meron or antimeron located at the core region of the quantized 
flux. 

Next let us consider the ac Josephson problem. The time variation of the phase of a wave 
function plays a role of the chemical potential. When a Josephson junction is constructed 
with the present system, the time derivative of the phase difference Ax across the junction 
is related to the chemical potential difference across the junction; namely, the voltage V 
across the junction is related to Ax as 

l -A X = qV/h. (52) 

The solution to the above equation is 

Ax = -j-t + const. (53) 

Since x an d X + 2vr are physically the same, the above relation implies the current across 
the junction will oscillate with the frequency 2qV/h. Overall, the flux quantization in the 
units y an d the Josephson frequency are both explainable in the present model. 

The phase of the superconducting order parameter is currently believed to stem from the 
breakdown of the global 17(1) gauge invariance^ and it is also believed the particle number 
nonconserving formalism is necessary to take into account it. However, the essential features 
of superconductivity are obtainable if a phase mode that parameterizes U{l)/Z 2 exists;— 
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thus, superconductivity without the breakdown of the global U(l) gauge invariance is in 
principle possible. Indeed, the present theory shows such a phase, |% arises spontaneously 
without the breakdown of the global U(l) gauge invaraince. 

VIII. CONCLUSION 

We have considered persistent current generation by small polarons embedded in a spin- 
density wave background in a two dimensional lattice. A spin current is generated by a spin 
Berry phase arising from a spiral spin texture produced by merons or antimerons. This state 
is currentless due to the exact cancellation of the currents for up- and down-spin electrons. 
However, when the polarons form side-by-side bipolarons, the motion of the system as a 
whole occurs, spontaneously; and the persistent current appears with the bipolaronic cores 
of merons and antimerons serving as pumps for it. The present model explains qualitatively 
various anomalous features of the cuprate superconductors. It also suggests the appearance 
of the phase variable or "Nambu-Goldstone" mode for the superconductivity without the 
breakdown of the global U(l) gauge invariance. 
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